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INTRODUCTION 


There  are  many  situations  where  it  is  useful  to  know  the  probability  that  some 
weather  event  will  occur  on  some  specified  date  in  the  future.  One  way  to  estimate 
this  probability  is  to  examine  past  records,  and  find  the  proportion  of  time  that 
the  event  occurred  in  previous  years.  For  example,  if  one  were  interested  in  the 
probability  of  rain  in  Honolulu  on  May  8,  1982,  one  method  would  be  to  examine  the 
available  records  of  Honolulu  for  May  8,  for  previous  years.  Presently,  large 
environmental  data  baser  exist  which  can  be  used  to  estimate  the  probabilities  of 
various  weather  occurrences.  Because  of  the  massive  amount  of  data,  summaries  are 
usally  utilized.  USAFETAC,  Air  Weather  Service  regularly  produces,  for  a large 
number  of  stations.  a "Revised  Uniform  Summary  of  Surface  Weather  Observations." 
These  "RUSSWO's"  provide  a very  convenient  summary  and  are  widely  utilized. 

There  are  a number  of  situations  where  even  these  summaries  are  unwieldly. 

For  example,  decision  models  are  now  widely  used  in  a variety  of  disciplines.  In 
a number  of  these^ environmental  conditions  are  required  as  inputs.  To  input  raw 
historical  data  for  a decision  model  into  a computer  may  require  a prohibitive 
amount  of  storage  space.  If  one  can  give  frequency  distributions  for  the  states 
of  a given  weather  element,  then  it  may  be  possible  to  substitute  models  (formulas) 
for  a summary.  This  would  substantially  increase  the  number  of  situations  where 
it  would  be  possible  to  utilize  decision  models  requiring  environmental  inputs. 

Part  I of  this  report  gives  statistical  models  for  the  probability  of  daily 
precipitation  amounts.  Part  II  gives  models  for  the  probability  of  "extreme 
values"  for  precipitation  during  a given  month.  In  all  cases  the  models  are  de- 
signed to  answer  questions  requiring  as  answers,  estimation  of  the  cumulative 
distribution  function.  This  report  gives  models  for  ten  stations  at  diverse  world- 
wide locations.  In  addition  to  separate  models  for  each  month,  additional  models 
are  given  for  which  month  of  the  year  is  an  input  to  a general  model,  valid  for 
all  months.  For  the  latter  model  if  one  wished  results  for  May,  for  example,  we 
would  set  M=5,  and  similar ily  M=ll  for  November.  As  would  be  expected,  the  general 
model  does  net  give  results  as  accurate  as  the  individual  models.  An  estimate  of 
the  accuracy  of  the  formulas  is  given  in  all  cases.  The  models  were  developed 
from  data  obtained  from  a "Revised  Uniform  Summary  of  Surface  Weather"  for  the 
appropriate  station,  published  by  USAFETAC,  Air  Weather  Service. 
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PART  I DAILY  PRECIPITATION 


l.  The  Lognormal  Model  For  Dally  Precipitation  I 

l 

Let  the  proportion  of  days  when  there  is  no  rain  be  denoted  by  D,  and  let  , 

F(x)  be  the  probability  of  an  amount  of  rain  less  than  or  equal  to  x.  Assume 

for  the  days  when  rain  jccurs,  that  the  amount  of  rain  has  a lognormal  dis-  , 

tribution.  Then,  if  4>(z)  is  the  cumulative  normal  distribution,  we  may  write 

F(x)  = D + (l-D)<J>(z)  . I 

where  z = (Zn  x-p)/o,  and  p and  o are  the  mean  and  standard  deviation  of  Zn  x.  | 

For  any  specified  month,  the  lognormal  model,  denoted  by  MODEL  I,  was  found 
to  fit  the  data  quite  well  for  each  station  for  each  month.  Two  different  methods 
were  used  to  estimate  the  parameters  p,  o and  D.  The  first  and  simplest  was  to  ^ 

estimate  D as  the  proportion  of  days  when  no  rain  occurred.  Estimates  of  p and  o ^ 

were  then  obtained  using  the  method  of  maximum  likelihood.  The  maximum  likelihood 
estimates  of  p and  a are  the  arithmetic  mean  and  the  standard  deviation  (with 
divisor  the  sample  size  n instead  of  n-1)  of  the  natural  logarithms  of  the  data 

for  days  when  there  was  rain.  i 

The  second  method  obtained  estimates  of  p,  o and  D using  least  squares  non- 
linear regression  techniques.  Let  y = C(x)  be  the  proportion  of  the  days  when 
the  rainfall  is  x inches  or  less.  Then  we  regress  y on  F(x)  = D + (1-D)4>(z).  The 
statistical  package  SAS  has  a non-linear  regression  procedure  called  NONLIN,  and 
this  was  used  to  obtain  estimates  m,  s and  d for  p,  o and  D.  This  technique 
minimizes  the  sum  of  squares  between  the  empirical  cumulative  distribution  function 
and  the  model  cumulative  distribution  function  at  the  boundary  points  of  the  x- 
intervals . 

Since  we  are  in  general  interested  in  estimating  the  probabilities  of  the 
occurrences  of  rainfall  less  than  (greater  than)  some  specified  amount  x,  and  not 
merely  the  estimates  of  p , o and  D,  the  second  procedure  seems  the  more  appropriate. 

Section  5 gives  the  estimates  m,  s and  d from  the  second  method  for  each  of  the 
stations  for  each  month  of  the  year. 

Model  1 results  in  separate  values  of  m,  s,  and  d for  each  month  of  the  year. 

Models  II  and  III  are  comprehensive  models  for  a given  station.  The  twelve 
monthly  values  for  each  of  m,  s and  d in  Model  II  are  represented  by  a cubic 
polynomial  in  M,  the  month  of  the  year  (M=l  for  January,  M=2  for  February,  etc.). 

Model  III  was  obtained  by  examining  the  residuals  from  Model  II  and  "fitting" 
the  residuals  to  an  appropriate  sine  term.  The  twelve  monthly  values  for  each  of 

m,  s and  d in  Model  III  are  thus  represented  by  a cubic  polynomial  in  M plus  a 
sine  term. 

The  polynomials  for  each  of  Models  IT  and  III  are  given  in  Section  5. 
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2 . The  Kappa  and  the  Modified  Kappa  Mod  el  For  Da i ly  Precipi tat ion 


Mielke  (1973)  has  recently  introduced  a family  of  frequency  distributions. 

This  family  has  been  termed  the  Kappa  distribution  and  its  cumulative  distribution 
function  is  given  by 

F(x)  = [(x/B)a0/(a  + (x/3)a°)  l1/(l°  a , 3 , 0 > 0. 

For  precipitation,  Mielke  used  0=1,  and  termed  the  resulting  family  the  "two 
parameter  Kappa  distribution".  For  the  data  from  the  stations  studied  in  this 
report,  both  the  "three  parameter  Kappa  model"  and  the  "modified  three  parameter 
model"  gave  fits  which  were  better  than  those  for  the  "two  parameter"  family. 

The  "modified  Kappa  model"  has  a cumulative  distribution  given  by 

F(x)  = [ (x/B)  lt)/ (a  + (x/B)a0)  11/,a  a,B,0  > 0. 

Essenwanger  (1976)  has  given  the  maximum  likelihood  estimates  for  a , 3 and  0 for 
the  Kappa  model.  However,  for  our  modeling  purposes  we  prefer  the  estimates 
obtained  by  the  "least  squares"  method.  Using  this  method,  the  estimates  n,b,  ; 
for  u,  3 and  0 are  those  values  which  minimize  the  sum  of  squares  between  the 
observed  empirical  cumulative  distribution  function  and  the  model  cumulative  dis- 
tribution function.  The  empirical  cumulative  distribution  function  was  obtained 
from  the  appropriate  RUSSWO's.  The  categories  of  "none"  and  "trace"  were  com- 
bined to  give  "no  precipitation".  "No  precipitation"  is  taken  to  mean  precipita- 
tion less  than  or  equal  to  .005  inches. 

Two  procedures  were  used  to  obtain  the  "least  squares  estimates".  One  was 
non  linear  regression  of  the  observed  empirical  cumulative  distribution  values 
on  the  model  cumulative  distribution  function.  The  SAS  (Statistical  Analysis 
System)  procedure  NONLIN  was  used.  An  equivalent  method  was  to  use  an  unsophisti- 
cated "grid  search".  Using  the  data  from  the  RUSSWO's,  the  sum  of  squares  of  the 
differences  between  the  empirical  cumulative  distribution  function  and  the  model 
cumulative  distribution  function  was  calculated  for  a grid  of  values  for  a,  3 and 
0.  The  "least  squares"  estimates  a,  b,  t are  those  values  of  a,  3 and  0 in  the 
grid  for  which  the  above  sum  of  squares  of  differences  was  smallest. 

In  Section  5,  we  give  the  values  a,  b,  t for  the  "Modified  Kappa"  model.  We 
term  this  Model  IV.  A comprehensive  model  for  each  station  was  also  developed 
using  the  "Modified  Kappa"  model.  In  Model  V,  the  values  a,  b,  t are  given  as 
quadratic  polynomials  in  M,  the  month  of  the  year.  The  polynomials  for  Model  V 
are  given  in  Section  5. 

It  should  be  noted  that  to  obtain  the  est imated  probabi lity  of  no  precipitation 
(including  trace)  using  the  modified  Kappa  model,  we  put  x =.005  inches. 

3 . Goodness  of  Fit  of  The_  Lognormal  and  Modified  Kappa  Models 

It  is  expected  that  the  models  will  be  used  for  estimating  the  probabilities 
of  rainfall  less  than  (greater  than)  a specified  amount.  Thus  the  measure  of 
"goodness"  used  was  the  root  mean  square  of  the  difference  between  the  empirical 
cumulative  distribution  function  and  the  model  cumulative  distribution  function, 
summed  over  all  the  boundary  points  between  categories  in  the  table  of  observed 
values  in  the  RUSSWO's. 
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As  previously  mentioned,  the  three  parameter  modified  Kappa  distribution  gave 
decidedly  better  fits  than  the  two  parameter  kappa  distribution.  The  three  para- 
meter modified  Kappa  gives,  on  the  average,  slightly  better  fits  than  the  lognormal 
distribution.  An  advantage  of  either  Kappa  distribution  is  that  no  tables  are  needed 
to  calculate  needed  probabilities  since  each  is  a closed  form  distribution. 

Section  5 gives  the  RMS  of  the  fits  for  all  of  the  lognormal  and  modified 
Kappa  distributions. 

4.  Use  of  the  Models 


Suppose  we  wish  to  estimate  the  probability  of  more  than  0.5  inch  of  rain  on 
June  20,  1984,  in  Bedford,  Massachussetts.  We  first  find  the  probability  of  0.5 
inch  or  less. 

Lognormal  model  for  June 

Putting  x = .5,  we  have 

z = (£n  x-y) /a 

= (-.693— 3.336)/2.230 
= 1.185 

then  the  probability  of  0.5  inch  or  less  is  estimated  by 
D + ( 1-D)  $ (z)  = .532  + (1-.532)  4>  (1.185) 

= .945  . 

The  probability  of  more  than  0.5  inches  is  1-.945  = .055. 

Kappa  model  for  June 

Putting  x = .5,  and  using  a = 70,  b = .450  and  t = .08 
F(x)  = [(.5/.450)5’6  / (70  + (.5/.450)5-6)]1/7° 

= .949 

The  probability  of  more  than  0.5  inches  is  estimated  to  be  .051. 

Lognormal  models  (overall  - Model  II) 

From  Section  5,  using  M = 6,  we  calculate 
d = .504,  m = -3.923,  s = 2.314. 

Then  z = (-.  693— 3 . 923) /2 . 314 
= 1.396. 

Estimated  probability  of  less  than  .5  inches  is  estimated  by 
.504  + (1-.504)  $ (1.396)  = .960. 

The  probability  of  more  than  .5  inches  is  estimated  to  be  .040. 


Tables  of  Models  for  Dai ly  Precipitation 
Station  12867  Patrick  AFB 
Model  I (Lognormal) 


d 

in 

s 

RMS  of  Fit 

.Jan. 

.663 

-3.757 

2.321 

.022 

Feb. 

.618 

-3.233 

2.370 

.024 

Mar. 

. 649 

-3.111 

2.326 

.018 

Apr . 

.714 

-3.445 

2.392 

.018 

May 

.611 

-3.130 

2.383 

.023 

Jun . 

. 494 

-2.620 

2.318 

.025 

Jul  . 

.527 

-2.845 

2.295 

.026 

Aug. 

.504 

-2.915 

2.304 

.022 

Sep. 

.384 

-2.675 

2.341 

.028 

Oct. 

. 466 

-3.046 

2.351 

.023 

Nov . 

.632 

-3.743 

2.116 

.018 

Dec . 

.656 

-3.954 

2.201 

.021 

Model  11  (Lognormal) 

d = .557  4 . 105M  - 

.027M2  + .002MJ 

m = -3.712  + 

. 07  5M  - 

.031M2  - .OOSM^ 

s = 2.314  + 

. 021M  - 

.002M2 

RMS  Error  = .128 

Model  III  (Lognormal) 


.638 

+ 

. 054M  - 

.019M2  + 

.001M3  - 

.050 

sin (n (M- . 5) /2) 

-3.950 

+ 

. 224M  + 

.005M2  - 

•002M}  + 

.145 

sin(n (M- . 5 ) / 2 ) 

2 . 309 

+ 

. 02 3M  - 

. 003M2 

- 

.026 

sin (n (M-L  5 ) / 2 ) 

RMS  Error 


Station  12867  Patrick  AFB  (CONT.) 


Model  IV  (Modified  Kappa) 


a 

b 



RMS  of  Fit 

Jan. 

40 

. 100 

.05 

.015 

Feb. 

no 

.975 

.05 

.011 

Mar . 

100 

.800 

.05 

.009 

Apr . 

30 

.100 

.05 

.017 

May 

110 

.975 

.05 

.010 

Jun. 

100 

1.325 

.08 

.013 

Jul. 

80 

.800 

.08 

.016 

Aug. 

110 

1.150 

.08 

.008 

Sep. 

40 

.975 

.11 

.010 

Oct . 

60 

.800 

.08 

.008 

Nov . 

70 

.275 

.05 

.006 

Dec . 

80 

.275 

.05 

.007 

Model  V (Modified  Kappa) 

a = 61.59  + 

7.880M  - 

.652M2 

b = .0095  + 

. 361M  - 

.029M2 

t = .021  + 

. 015M  - 

.001M2 

RMS  Error  = 0.036 
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Station  13802  Scott  AFB 


Model  1 (Lognormal) 


d 

m 

s 

RMS  ol  Hi t 

Jan . 

.520 

-4.016 

2.198 

.031 

Feb. 

.489 

-3.825 

2.275 

. 004 

Mar. 

.466 

-3.504 

2.272 

.003 

Apr . 

.467 

-3.078 

2.341 

.036 

May 

.517 

-2.970 

2.387 

.029 

Jun . 

. 562 

-2.924 

2.370 

.024 

JuJ  . 

.613 

-2.994 

2.377 

.023 

Aug. 

. 663 

-3.120 

2.490 

.021 

Sep . 

.658 

-3.084 

2.403 

.020 

Or  t . 

.667 

-3.080 

2.273 

.017 

Nov . 

.602 

-3.344 

2.303 

.025 

Dec . 

. 526 

-3.661 

2.295 

.003 

Mode!  11  (Lognormal) 

d = .610  - , 

. 1 I6M  + 

.027M2  - .002M3 

m = -4.465  + , 

,444M  - 

.035M2 

s = 2.113  + 

. 074M  - 

. 004M2 

RMS  Error  = .121 

Model  III  (Lognormal) 

d = .610  - . 

, 1 I6M  + 

.027M2  - .002M3 

m = -4.519  + . 

. 49 1 M - 

.043M2  + .001M1  - 

.093  sin 

( tt  (M—  1. 5)/2-5) 

s = 2.108  + . 

. 075M  - 

. 005M2 

.027  sin 

(n (M-l. 5) /2) 

RMS  Error  = .030 
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Station  13802  Scott  AFB  (CONT.) 


Model  IV  (Modified  Kappa) 


a 

b 

t 

RMS  of  Fit 

Jan. 

110 

.625 

.05 

.010 

Feb. 

100 

.975 

.05 

.018 

Mar. 

80 

.625 

.08 

.011 

Apr . 

100 

1.150 

.08 

.021 

May 

120 

.975 

.08 

.016 

Jun. 

40 

.450 

.08 

.018 

Jul . 

120 

1.325 

.05 

.015 

Aug. 

90 

.625 

.05 

.013 

Sep. 

120 

.975 

.05 

.010 

Oct. 

90 

.625 

.05 

.009 

Nov. 

120 

1.150 

.05 

.013 

Dec . 

110 

.975 

.05 

.014 

Model 

V (Modi 

fied  Kappa) 

a = 

110. 

- 6.409M  + 

.584M2 

b = 

.756 

+ . 023M  - 

•001M2 

t = 

.053 

+ . 006M  - 

.001M2 

RMS  Error  = 0.024 
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Station  14601  Bangor,  Maine 


Model  I (Lognormal) 


d 

m 

S 

RMS  of  Fit 

Jan. 

.436 

-3.135 

2.259 

.037 

Feb. 

.436 

-3.158 

2.265 

.034 

Mar . 

.448 

-3.274 

2.227 

.032 

Apr. 

.443 

-3.312 

2.158 

.031 

May 

.450 

-3.349 

2.196 

.030 

Jun. 

.463 

-3.370 

2.269 

.032 

Jul . 

.476 

-3.390 

2.723 

.028 

Aug. 

.516 

-3.473 

2.172 

.024 

Sep. 

.530 

-3.203 

2.295 

.026 

Oct . 

.520 

-3.060 

2.272 

.022 

Nov. 

.416 

-2.969 

2.399 

.040 

Dec. 

.425 

-3.353 

2.224 

.027 

Model  II  (Lognormal) 

d = .487  - 

. 052M  + 

.013M2  - .001M3 

m = -2.834  - 

. 238M  + 

.032M2  - .001M3 

s = 2.354  - 

. 09  3M  + 

.016M2  - .001M3 

RMS  Error  = .155 

Model  III  (Lognormal) 

d = .476  - 

. 04 5M  + 

•011M2  - .001M3  + 

.020  Sin 

(u (M-l . 5)  / 3) 

m = 2.873  - 

. 210M  + 

.026M2  - .001M3  + 

,073  Sin 

(n(M-1.5)/3) 

s = 2.364  - 

. 097M  + 

.017M2  - .001M3  + 

.049  Sin 

(it  (M-l . 5)  / 2) 

* j RMS  Error  = .034 
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i 


I 

I 


Station  14601  Bangor,  Maine  (CONT.) 
Model  IV  (Modified  Kappa) 


a 

b 

t 

RMS  of  Fit 

Jan. 

40 

.450 

.11 

.024 

Feb. 

40 

.450 

. 11 

.022 

Mar , 

110 

.975 

.08 

.017 

Api . 

50 

.450 

.11 

.015 

May 

90 

.800 

.08 

.016 

Jun. 

100 

.800 

.08 

.013 

Jul . 

50 

.450 

.08 

.007 

Aug. 

40 

.275 

.08 

.009 

Sep. 

90 

.625 

.08 

.014 

Oct . 

40 

.450 

.08 

.010 

Nov . 

90 

1.675 

.08 

.025 

Dec . 

110 

.975 

.08 

.014 

Model 

y_  (Mod if 

led  Kappa) 

a = 

56.136  + 

1.741M 

+ . 

062M2 

b = 

.724  - 

. 0S7M 

+ . 

01 0M2 

t = 

.119  - 

. 009M 

+ . 

0005M2 

RMS  Error  = 0.020 


1-9 


Station  14702  Bedford,  Mass ac h u s s . ts 


Model  I (Lognormal) 


d 

m 

S 

RMS  of  Fit 

Jan. 

.467 

-3.097 

2.322 

.031 

Feb . 

.450 

-3.214 

2.271 

.030 

Mar. 

.487 

-3.129 

2.337 

.032 

Apr. 

.454 

-3.052 

2.271 

.032 

May 

.487 

-3.159 

2.233 

.029 

Jun . 

.532 

-3.336 

2.230 

.026 

Jul. 

.542 

-3.461 

2.308 

.027 

Aug. 

.570 

-3.301 

2.354 

.025 

Sep. 

.591 

-3.305 

2.323 

.020 

Oct. 

.604 

-3.057 

2.320 

.019 

Nov . 

.492 

-3.104 

2.469 

.032 

Dec  . 

.488 

-3.241 

2.371 

.030 

Model 

1 d = 

11  (Lognormal) 

.316  - . 062M  + 

.016M2  - .001M3 

l 

-3.101 

- . 034M  + 

.001M2 

1 

s = 

2.356 

- . 04 3M  + 

•006M2 

RMS  Error  = .140 

Mode  1 

d = 

III  (Lognormal) 

.531  - . 072M  + 

.018M2  - .001M3  - 

.011  Sin 

(it  (M-  . 5)  /1 . 5) 

m = 

-2.956 

- . 121M  + 

.016M2  - .001M3  - 

.077  Sin 

(if(M-  .5)/3) 

1 s = 

2.431 

- . 090M  + 

. 014M2  - . 001M3  + 

.046  Sin 

(tt  (M-2 . 5)  /2) 

RMS  Error  = .030 
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Station  14702  Bedford,  Massachussets  (CONT , ) 


Model  IV  (Modified  Kappa) 


a 

b 

t 

RMS  of  Fit 

Jan. 

120 

1.150 

.08 

.017 

Feb. 

90 

.975 

.08 

.014 

Mar. 

100 

.975 

.08 

.014 

Apr. 

120 

1.325 

.08 

.021 

May 

90 

.975 

.08 

.015 

Jun. 

70 

.450 

.08 

.013 

Jul. 

120 

1.325 

.05 

.015 

Aug. 

30 

.275 

.08 

.018 

Sep. 

110 

.975 

.05 

.011 

Oct. 

120 

1.325 

.05 

.013 

Nov, 

90 

.975 

.08 

.017 

Dec . 

90 

.800 

.08 

.013 

Model  V (Modified  Kappa) 

a = 129.0  - 13.045M  + .894M2 
b = 1.327  - . 145M  + .010M2 

t = .093  - . 006M  + .0004M2 

RMS  Error  = 0.020 


^ 


t 

I 


Station  22521  Honolulu,  Hawaii 
Model  I (Lognormal) 


d 

m 

s 

RMS  of  Fit 

Jan. 

.448 

-3.830 

2.315 

.034 

Feb. 

.410 

-4.097 

2.103 

.032 

Mar. 

.385 

-4.385 

2.094 

.044 

Apr. 

.342 

-4.729 

1.837 

.044 

May 

.385 

-5.034 

1.637 

.047 

Jun. 

.347 

-5.267 

1.279 

.042 

Jul . 

.355 

-5.194 

1.342 

.042 

Aug. 

.340 

-5.070 

1.491 

.005 

Sep. 

.411 

-4.934 

1.590 

.040 

Oct . 

.376 

-4.649 

1.871 

.032 

Nov . 

.344 

-4.479 

2.021 

.045 

Dec . 

.320 

-4.234 

2.139 

.046 

Model  II  (Lognormal) 

d = .525 

- . 086M 

+ 

.013M2  - .001M3 

m = -3.132 

- . 6 1 3M 

+ 

.052M2  - .001M3 

s = 2.674 

- . 303M 

+ 

.014M2  + .001M3 

RMS  Error  = .105 

Model  III  (Lognormal) 

d = .534 

- . 090M 

+ 

.014M2  - .001M3  - 

.012  Sin 

(it  M/1 . 75) 

m = -3.203 

- . 562M 

+ 

.041M2  + 

.134  Sin 

(n(M-l. 5)/3) 

s = 2.842 

- . 394M 

+ 

.028M2  + 

.124  Sin 

(ir  (M-2 . 5)  /3) 

RMS  Error  = .044 
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Station  22521  Honolulu,  Hawaii  (CONT. ) 


Model  IV  (Modified  Kappa) 


a 

b 

t 

RMS  of  Fit 

Jan. 

20 

.100 

.08 

.012 

Feb. 

20 

. 100 

co 

o 

.005 

Mar. 

110 

.625 

.05 

.014 

Apr 

40 

.100 

.08 

.010 

May 

70 

.100 

.05 

.007 

Jun. 

120 

.100 

.05 

.012 

Jul. 

80 

.100 

.05 

.011 

Aug. 

70 

.100 

.05 

.011 

Sep. 

70 

.100 

.05 

.010 

Oct. 

30 

. 100 

.08 

.011 

Nov. 

30 

.100 

.08 

.014 

Dec. 

20 

. 100 

.08 

.005 

Model  V (Modified  Kappa) 

a = -5.909  + 29.066M  - 2.333M2 

b = .100 

t = .096  - . 014M  + .001M2 

RMS  Error  = 0.035 
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Si-ition  26435  Nenana,  Alaska 


Model  I (Lognormal) 


d 

m 

s 

RMS  of  Fit 

Jan. 

. 565 

-4.317 

1.649 

.018 

Feb. 

.537 

-4.188 

1.610 

.025 

Mar . 

.611 

-4.450 

1.572 

.020 

Apr . 

.653 

-4.557 

1.552 

.018 

May 

.633 

-4.325 

1.807 

.021 

Jun. 

.477 

-3.687 

2.004 

.028 

Jul . 

.428 

-3.391 

1.976 

.029 

Aug. 

.390 

-3.446 

1.988 

.025 

Sep . 

.498 

-3.709 

1.863 

.022 

Oct . 

.491 

-4.275 

1.624 

.024 

Nov . 

.506 

-4.330 

1.689 

.022 

Dec . 

.486 

-4.613 

1.460 

.024 

Model  II  (Lognormal) 

d = .478  + . 095M 

- .021M2  + .001M3 

m = -4.011 

- . 388M 

+ . 109M2  - . 007M  3 

s = 1.653 

- . 08 5M 

+ . 034M2  - .002M3 

RMS  Error  = .126 

Model  III  (Lognormal) 

d = .578 

+ . 027M 

- .009M2 

.065  Sin 

(n  (M- . 5)  / 3) 

m = -4.341 

- . 166M 

+ .066M2  - .004M3  + 

.254  Sin 

(n(M-.5)/3) 

s = 1.705 

- . 08 5M 

+ .031M2  - .002M3  - 

.038  Sin 

(n (M-l)/3) 

RMS  Error  = .034 
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Station  26435  Nenana,  Alaska  (CQNT. ) 


Model  IV  (Modified  Kappa) 


a 

b 

t 

RMS  of  Fit 

Jan . 

50 

. 100 

.08 

.009 

Feb. 

}0 

. 100 

.08 

.009 

Mar. 

120 

. 100 

O 

00 

.011 

Api . 

120 

. 100 

.05 

.010 

May 

60 

.100 

.05 

.004 

Jun . 

100 

.450 

.08 

.012 

Jul. 

90 

.450 

.11 

.015 

Aug. 

70 

.450 

. 1 1 

.012 

Sep. 

30 

.100 

.11 

.015 

Oct. 

70 

.100 

. 11 

.009 

Nov. 

110 

.100 

.11 

.010 

Dec. 

110 

. 100 

.11 

.009 

Model  V (Modified  Kappa) 

a = 72.27  + 1.404M  + .005M2 
b = - .099  + . 116M  - .0087M2 
t = .070 

RMS  Error  = 0.045 
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Station  33123  Tripoli,  Libya 


Model  I (Lognormal) 


d 

m 

s 

RMS  of  Fit 

Jan. 

.565 

-3.492 

2.150 

.027 

Feb. 

.704 

-3.793 

2.086 

.017 

Mar. 

. 782 

-4.136 

2.109 

.016 

Apr. 

.794 

-4.623 

1.910 

.020 

May 

.842 

-5.051 

1.614 

.013 

Jun. 

.900 

-5.517 

1.076 

.007 

Jul. 

.975 

-5.358 

1.128 

.002 

Aug. 

97? 

-4.956 

1.495 

.002 

Sep. 

.853 

-4.627 

1.957 

.017 

Oct . 

.630 

-4.045 

2.129 

.025 

Nov . 

.590 

-3.607 

2.290 

.026 

Dec . 

.545 

-3.040 

2.253 

.026 

Model  ii  (Lognormal) 

G. 

II 

o 

+ . 088M 

- .001M2  - .001M3 

m = -3.024 

- . 565M 

+ .023M2  + . 002M3 

s = 2.452 

- . 257M 

+ .OHM2  + .001M3 

RMS  Error  = .074 

Model  III  (Lognormal) 

d = .301  + 

. 204M 

-.019M2  +. 

070  Sin 

(n(M-  . 5) /3) 

m = --2.876 

- . 531M 

+.003M2  + .004M3  +. 

213  Sin 

(it  (M-l . 5)  / 3) 

s = 2.810 

- . 425M 

+.031M2  +. 

334  Sin 

(it  (M-2  )/3) 

RMS  Error  = .024 
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Model  IV  (Modified  Kappa) 


f 


(CONT.) 


a 

b 

t 

RMS  of  Fit 

Jan. 

60 

.275 

.08 

.018 

Feb. 

120 

.275 

.05 

.009 

Mar. 

120 

.275 

.02 

.014 

Apr . 

110 

. 100 

.02 

.007 

May 

120 

. 100 

.02 

.023 

Jun. 

120 

. 100 

.02 

.036 

Jul. 

120 

.100 

.02 

.040 

Aug. 

120 

.100 

.02 

.038  ; 

Sep. 

120 

.100 

.02 

.020 

Oct. 

40 

. 100 

.05 

.009 

Nov. 

100 

.625 

.05 

.010  : 

Dec. 

120 

.800 

.08 

.015 

* 

Model  V 

a = 78 

(Modified  Kappa) 

.409  + 11.941M  - 

.927M2 

b = .569  - . 187M  + .0164M2 

t = .094  - . 026M  + , 002M2 

RMS  Error  = 0.022 
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Station  41108  Saigon,  Vietnam 


Model  I (Lognormal) 


d 

m 

s 

RMS  or  Fir 

Jan. 

.859 

-4.329 

1.907 

.009 

Feb. 

.961 

-4.371 

1.836 

.003 

Mar. 

.879 

-4.104 

1.837 

.005 

Apr . 

.714 

-4.061 

2.033 

.O1^ 

May 

.317 

-2.060 

2.002 

.034 

Jun. 

. 167 

-2.027 

2.134 

.047 

Jul. 

.140 

-2.010 

2.123 

.044 

Aug. 

. 129 

-1.890 

2.041 

.044 

Sep. 

.100 

-1.924 

2.078 

.040 

Oct . 

.207 

-2.021 

2.066 

.033 

Nov. 

. 500 

-3.066 

2.248 

.026 

Dec . 

.755 

-3.428 

2.129 

.014 

Model  ll  (Lognormal) 

d = .889 

+ . 099M 

- .058M2  + .004M3 

m = -4.667 

+ . 019M 

+ .117M2  - .009M3 

s = 1.813 

+ . 03 1M 

+ .003M2 

RMS  Error  = .224 

Model  III  (Lognormal) 

d = .834  + .139M 

- .067M2  + .005M3  + 

.103  Sin 

(ir(M-1.5)/3) 

m = -4.511 

- . 094M 

+ .141M2  - .OHM3  - 

.291  Sin 

(it  (M-l . 5)  / 3 ) 

s = 1.858 

- . 00 2M 

+ .010M2  - .001M3  - 

.084  Sin 

(ir(M-1.5)/3) 

RMS  Error  = .050 
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Station  41103  Saigon,  Vietnam  (CONT.) 
Model  IV  (Modif  ied  Kappa) 


a 

b 

t 

RMS  of  Fit 

Jan. 

120 

.100 

.02 

.007 

Feb . 

120 

. 100 

.02 

.033 

Mar. 

120 

.100 

.02 

.013 

Apr. 

80 

. 100 

.05 

.011 

May 

90 

1.325 

. 17 

.025 

Jun . 

30 

1.150 

.23 

.027 

Jul . 

50 

1.500 

.23 

.025 

Aug. 

50 

1.675 

.23 

.035 

Sep. 

40 

1.675 

.23 

.035 

Oct. 

20 

.800 

.23 

.020 

Nov. 

120 

.975 

.08 

.010 

Dec . 

50 

.100 

.05 

.016 

Model  V (Mod if 

ied  Kappa) 

a = 164.318 

- 27.585M  + 

1.646M2 

b = -.542 

+ . 607M  - 

.045M2 

t = -.079 

+ . 091M  - 

.0066M2 

RMS  Error  = 0.116 
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S t_a  lion  457 15  Shemya,  A 1 a ska 
Model  I (Lognormal) 


d 

m 

s 

RMS  of  Fit 

Jan 

.028 

-3.716 

1.770 

.035 

Feb. 

.017 

-3.986 

1.731 

.034 

Mar 

.050 

-3.992 

1.806 

,041 

Apr 

. 100 

-4.217 

1.820 

.005 

May 

. 124 

-4.217 

1.895 

.043 

Jun 

.114 

-4.491 

1.786 

.053 

Jul 

. 105 

-4.161 

2.031 

.055 

Aug 

. 147 

-4.018 

1.966 

.042 

Sep 

. 181 

-3.726 

1 .937 

.003 

Oet 

.076 

-3.426 

1.890 

.027 

Nov 

.042 

-3.449 

1.849 

.030 

Dec 

. 059 

-3.728 

1.879 

.043 

Mode  l 

11  (Lognormal) 

d = 

- .004  + . 01 9M 

+ 

002M2 

m = 

-3.007  - . 654M 

+ 

101M2  - ,004MJ 

s = 

1.726  - . 006M 

-f 

007M2  - . 00 1M  3 

RMS  Error  = .240 

Model 

111  (Lognormal) 

d = 

- .010  + . 02 1M 

+ 

002M2 

.028  Sin 

(ir(M-1.5)/2) 

m = 

-2.746  - . 795M 

+ 

122M2  - .005M5  - 

.192  Sin 

(ir(M+  . 5)  / 3) 

s = 

1.718  + .010M 

+ 

007M2  - .001M3  - 

.040  Sin 

(n(M-l . 5 ) / 2 ) 

RMS  Error  = .045 
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Station  457  1 5__Sliernya  , A 1 a ska  (CONT  . ) 
Model  IV  (Modified  Kappa) 


a 

b 

t 

RMS  of  Fit 

Jan. 

60 

.27  5 

.23 

.034 

Feb. 

120 

.275 

.23 

.028 

Mar. 

1 10 

.275 

.20 

.01  7 

Apr. 

20 

.100 

.20 

.017 

May 

50 

.275 

. 14 

.015 

dun . 

20 

. 100 

. 14 

.008 

Jul  . 

40 

.275 

.14 

.01 1 

Aug. 

30 

.275 

.14 

.014 

Sep . 

40 

.275 

.17 

.016 

Oct. 

20 

.275 

.23 

.030 

Nov . 

80 

.450 

.23 

.041 

Dec. 

50 

.275 

.23 

.013 

Model  V (Modi 

a = 120. 455 

b = .303 

f ied  Kappa) 

- 22.273M  + 

. 03 5M  + 

.434M': 

.0034M2 

t = .287 

. 042M  + 

. 0033M ^ 

KMS  Error  = 0,033 


Station  93780  Christchurch,  New  Zealand 


Model  I (Lognormal) 


d 

m 

s 

RMS  of  Fit 

Jan. 

.559 

-3.471 

1.984 

.016 

Feb. 

.591 

-3.460 

1.972 

.016 

Mar . 

.571 

-3.175 

1.926 

.002 

Apr. 

.573 

-3.145 

1.961 

.014 

May 

.527 

-3.225 

1.998 

.018 

Jun . 

. 587 

-3. 368 

1.970 

.001 

Jul. 

.500 

-3.165 

1.940 

.002 

Aug. 

.595 

-3.261 

1.995 

.015 

Sep. 

.567 

-3.407 

1.814 

.015 

Oct. 

.616 

-3.699 

1.877 

.018 

Nov. 

. 590 

-2.910 

1.981 

.014 

Dec . 

. 624 

-3.355 

2.044 

.018 

Model  II  (Lognormal) 

d = .577 

- .G10M  + 

.001M2 

m = -3.727 

+ . 246M  - 

•037M2  + .002M3 

s = 1.885 

+ . 079M  - 

•017M2  + .001M3 

RMS  Error  = .132 

Model  III  (Lognormal) 

1 

d = .606 

- . 028M  + 

.004M2  + 

.023  Sin 

(it  (M-l . 5)  ) 

m = -3.631 

+ . 184M  - 

.025M2  + .001M3  - 

.094  Sin 

(u (M-l  )/2) 

s = 1.912 

+ . 059M  - 

.013M2  + .001M3  - 

.051  Sin 

(ir  (M-l . 5)  /3) 

RMS  Error  = .019 
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Station  93780  Christchurch,  New  Zealand  (CONT. ) 


Mock' I IV  (Mod  if icd  Kappa) 


a 

b 

t 

RMS  ol 

dan. 

60 

.276 

.08 

.00  7 

Feb, 

80 

.276 

.08 

. 009 

Mar. 

; oo 

.460 

.08 

. 009 

Apr. 

80 

. 4 60 

.08 

.010 

May 

60 

.450 

.08 

.013 

Jun . 

70 

.276 

.08 

.010 

dul . 

40 

.275 

. 1 1 

.012 

Aug. 

60 

.275 

.08 

.008 

Sep. 

60 

.275 

.08 

.011 

Oct. 

40 

. 100 

.08 

.012 

Nov . 

70 

.450 

.08 

.008 

Dec. 

120 

.625 

.05 

.013 

Mode  1 

a = 

V (Mod i 1 i ed 

100. k55  - 13. 

Kappa) 

. 177M  + 1 .01  DM7 

b = . 430  - . 04 5M  + .0039M'’ 

t = . 080 

RMS  Error  = .021 
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PART  II  "EXTREME  VALUES"  OF  DAILY  PRECIPITATION 


1 . Models  of  "Extreme  Values"  for  Prec ipi ta t i on 

Two  sets  of  models  were  developed.  The  First  set  consists  of  twelve  models, 
one  for  each  month  of  the  year.  The  models  are  of  the  form 

G(x)  = exp  ( exp(-( (x-p) io) ) ) 

where  p and  a are  constants  with  o > 0.  C(x)  is  the  probability  that  for  the 
month  in  question,  the  maximum  (extreme)  daily  rainfall  will  be  less  than  or  equal 
to  x.  Section  2 gives  the  rationale  for  the  selection  of  this  model  and  discusses 
the  methods  used  to  get  the  estimates  m and  s for  p and  a respectively.  Section  6 
gives  the  tabled  values  of  m and  s for  the  stations  involved,  and  also  an  estimate 
of  the  accurary  of  the  model. 

The  second  set  of  models  are  val id  for  any  month  of  the  year.  The  same  form  is 
used  as  for  the  models  for  the  individual  months.  However,  m and  s do  not  appear 
explicitly.  Instead  m and  s are  each  either  a polynomial  in  M,  or  a polynomial  in 
M plus  a sine  term,  where  M represents  the  month  of  interest.  For  May  M = 5,  while 
for  September  M = 9.  Section  3 gives  further  details  on  the  derivation  of  the 
models . 


2.  Rationale  for  the  Models  Used 


Let  X be  the  amount  of  rair.  in  a day  selected  at  random  from  a specified  month. 

Then  X will  be  a random  variable  with  some  probability  density  function  f (x) . In 

general,  if  we  have  n observations  X1  , X^,  ...  , X on  a random  variable  X,  and 

X,..  is  the  ith  smallest  of  the  n observations,  then  the  cumulative  distribution 

function  of  X,  . is  given  by 
(n) 

G(x)  = [F(x)Jn 


where  F(x)  is  the  cumulative  distribution  function  of  X.  It  is  frequently  the  case 
that  F(x)  is  not  known.  However,  if  n is  large,  for  functions  f(x)  where  X is 
unbounded  and  F(x)  decreases  at  least  as  rapidly  as  the  exponential  function,  then 
G(x)  has  the  limiting  form 

G(x)  = exp  (-exp(-((x-p)/o))) 


where  p and  0 are  constants  with  a > 0.  This  distribution  is  sometimes  called  the 
"extreme  value"  distribution. 


The  parameters  p and  o may  be  estimated  from  data  by  several  methods.  It  can 
be  shown  that  the  maximum  likelihood  estimate  of  o is  given  by  the  solution  of  the 
following  equation: 


a 


x exp(-x . 


exp(-x  / 


0. 


n-i 


I 


A 

The  equation  may  be  solved  for  o by  an  iterative  procedure, 
estimate  of  fi  is  then  obtained  from 


P 


I 2 

n 1=1 


exp( 


-Xi«7] 


The  maximum  likelihood 


Another  method  of  obtaining  values  for  p and  a is  by  the  direct  regression  t>: 
the  empirical  cumulative  distribution  function  on  the  model  cumulative  distributicn. 


Both  methods  were  used  in  the  study.  The  second  method  resulted 
fits.  Instead  of  the  empirical  cumulative  distribution,  the  expected 
cumulative  distribution  for  the  i t_h  order  statistic,  i/(n+l)  was  used 
pendent  variable,  and  non-linear  regression  techniques  were  used.  The 
U and  o chosen  were  those  which  minimized 

i / (n+1 ) - oTI 


n 

i£i|_i' 


“(O' 


in  super ioi 
valu-  of  c!. 
as  t he  de- 
values for 


where  F(x^;p,o)  = exp(-?xp(-((x^-p)/a)))  . 
3 . Spec i f i cat  ion  of  the  Models 


Both  sets  of  models  use  the  form 


c(x)  = exp(-exp(-((x-u)/o))) 

where  p and  o are  constants,  and  C(x)  is  the  probability  that  the  maximum  (extreme) 
daily  rainfall  for  a specified  month  will  be  less  than  or  equal  to  x.  We  use  m 
and  s for  the  estimates  of  p and  a respectively.  For  the  first  set  of  12  models 
(one  for  each  month)  the  values  for  m and  s are  given  in  Section  7.  We  shall  call 
this  MODEL  I. 


The  second  set  of  models  are  each  valid  for  an  arbitrary  month.  Four  of  these 
models  were  developed.  We  shall  call  them  Models  II,  III,  TV  and  V.  In  each  case, 
both  m and  s are  functions  of  the  month  M.  For  Model  II,  m and  s are  each  a 
quadratic  polynomial  in  M,  while  for  Model  III  they  are  cubics.  Models  IV  and 
V are  quadratics  and  cubics  respectively  with  an  additional  sine  term. 


For  model  II,  we  have 


m 

o 

E 

11 

+ m M + 

m2M2 

s 

so 

+ SjM  + 

s2M2 

For  model 

hi 

, we  have 

m 

m0 

+ m^M  + 

m2M2 

+ 

m3M3 

s 

= so 

+ s M + 

s2M2 

+ 

s3M3  . 

Values  for  m^, 

m ^ , 

m2’  S0’ 

V ! 

S2 

for  Model  II  are  given  under  Model  II  in  Section 

I 


TT.-2 


» 


is 


-v-'r 


. b 


J 


6.  Values  for  m^,  m^,  m^ , m.^,  s^,  Sj  , s^,  s.^  for  Model  111  are  given  under  Model 

III  in  the  same  section.  The  tabulated  values  were  obtained  from  regression  an- 
alyses . 


Models  IV  and  V were  obtained  by  examining  the  residuals  from  Models  II  and 
III  and  adding  an  appropriate  sine  term  to  the  expressions  for  m and  s.  New  re- 
gressions were  then  run.  For  model  IV  we  have 


m = m_  + m,M  •+•  ">  M + .a  sin  c 
0 12  J m 


s = s + s,M  + s„M  + s-sin  c 
0 12  3 s 


For  model  V,  we  have 


m = m_  + m,M  + m„M  + m..M  + m.sin  c 
0 1 2 3 4 m 

2 3 

s = s„  + s.M  + s„M  + s„M  + s.sin  c 
0 1 2 3 4 s 

Values  for  the  m.'s,  s.'s,  c and  c are  given  in  Section  6 under  Models  IV  and  V. 
The  arguments  of  the  sm  functions  were  obtained  from  an  inspection  of  the  re- 
siduals . 


4 . Goodness  of  Fit  of  the  Models 

The  measure  of  "goodness"  used  was  the  root  mean  square  of  the  residuals  from 
the  regressions.  These  measures  are  Cabled  in  Section  6.  Choice  of  a model  should 
not  be  made  without  first  referring  to  the  table.  As  may  be  observed,  fits  for 
models  II,  III,  IV  and  V were  inferior  to  the  fits  of  Model  1 and  in  many  cases 
substantially  inferior.  This  should  be  expected  since  one  expression  is  used  for 
all  12  months  while  Model  I has  a separate  expression  for  each  month.  Models 
III  and  V give  better  fits  than  Models  II  and  IV  respectively,  and  Models  IV  and 
V give  better  fits  than  Models  II  and  III  respectively.  In  general,  the  models 
with  the  most  terms  (or  expressions)  give  the  best  fits. 

The  advantage,  of  course,  of  Models  II,  III,  IV  and  V is  that  in  each  case 
a single  expression  may  be  used  for  all  months. 

It  should  be  noted  that  in  many  cases  the  "extreme  value"  model  does  not 
provide  as  good  a fit  to  the  data  as  is  desired.  Thus,  more  work  is  needed  on 
models  to  replace  the  "extreme  value"  model. 

5.  Use  of  the  Models 


The  following  is  an  example  of  the  use  of  the  models. 

Problem:  What  is  the  probability  that  the  maximum  dai 1 y amount  of  rainfall  in 

Honolulu  in  February  will  be  less  than  2 inches? 

Solution : The  model  is 

G(x)  = exp(-exp(-((x-p)/o))) 

Using  Model  I,  our  estimates  of  p and  o are  m = .580  and  s = .662. 


I 

I 
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f 


I 

I 

Substituting  these  values  in  G(x)  , we  obtain  a probability  of  .89. 

Using  Model  II,  in  = 1.202  - . 331M  + . 025M^ 
s = 1.772  - . 463M  = .032M2 

where  M = 2.  Substituting,  we  have  m = .640,  s = .974  and  the  estimated  prob- 
ability is  .78. 


t 

t 


6.  Tables  of  Models  for  "Extreme  Values"  for  Prec Ipi tat i on 
Station  12867  Patrick  AF B 
Model  I 


m 

s 

___  RMS 

Jan , 

0.555 

0.833 

0,034 

Feb. 

0.913 

0.770 

0.051 

Mar 

0.  780 

0.757 

0.030 

Apr. 

0.584 

0.563 

0,030 

May 

0.885 

0.787 

0.025 

Jun . 

1.239 

0.964 

0,024 

Jul . 

1.018 

0.404 

0,033 

Aug. 

1.276 

0,  562 

0.039 

Sep. 

1.828 

1.034 

0.038 

Oct . 

1.201 

1 .155 

0.033 

Nov. 

0.653 

0.858 

0,048 

Dec . 

0.524 

0.479 

0.019 

Model 

ii 

m = 

.162  + 

. 293M  - .0206M2 

s = 

.750  + 

•0008M-  .0002M2 

RMS  Error 

= . 121 

Model 

III 

m = 

.975  - 

. 335M  + .0956M2  - 

,006M3 

s = 

1.173  - 

. 326M  + .0605M2  - 

.0031M* 

RMS  Error 

= .099 

Model 

JV 

. 309M  - .021M2  - . 

m = 

.057  + 

270 

Sin(Tr 

(M-2 . 5) / 2 ) 

s = 

.653  + 

. 016M  + .000M2  - . 

251 

Sin(ir 

(M-2 . 5 ) / 2 ) 

RMS  Error  = .092 

Model 

JV 

m = 

1.337  - 

, 565M  + .136M2  - 

.008M3  - 

,256  Sin(iT 

(M-3 , 5) /l . 5) 

s = 

1.194  - 

. 339M  + .063M2  - 

. 003M3  + 

.027  Sin(Tt 

(M-1.5)/l,5) 

RMS  Error  = .075 
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Station  13802  Scott  AKB 


Model  I 


m 

s 

RMS 

Jan. 

0.517 

0.462 

0.050 

Feb . 

0.613 

0.498 

0.026 

Mar. 

0.  752 

0.421 

0.023 

Apr. 

0.840 

0.517 

0.043 

May 

0.955 

0.591 

0.032 

Jun. 

1.095 

0.714 

0.036 

Jul. 

0.897 

0.734 

0.022 

Aug. 

0.938 

0.806 

0.027 

Sep. 

0.945 

0.  754 

0.031 

Oct . 

0.697 

0,569 

0.015 

Nov . 

0.722 

0.462 

0.028 

Dec . 

0.663 

0.603 

0.040 

Mode  l l l 

m = .328 

+ 

. 188M  - .014M2 

s = .279 

+ 

. 109M  - .007M2 

RMS  Error  = .050 

Model  1 I I 

tn  = .244 

+ 

. 253M  - .026M2  + 

.001M5 

s = .430 

— 

. 00 7M  + .014M2  - 

.001M3 

RMS  Error  = .049 

Model  VJ 

m = .328  + . 1 88M  - .014M2  + .001  Sln(it(M-1.5)) 

s = .214  + . 1 30M  - .009M2  - .095  Sin(tr  (M-2 . 5) /3) 

RMS  Error  = .046 

Model  V 

m = .239  + . 256M  - .026M2  + .001M3  - .004  Sin(tr  (M-l . 5) ) 

s = .310  + . 058M  - . 005M2  - .001M3  - .088  Sin(it(M-2. 5/3) 

RMS  Error  = .045 

II-  6 
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Station  14601  Bangor,  Ma i ne 


I 


Model  I 


m 

s 

KMS 

Jan. 

0.81  1 

0.371 

0.031 

Feb . 

».'»16 

0.285 

0.049 

Mar . 

0.  762 

0.  350 

0.034 

Apr . 

0.772 

0.2  58 

0.036 

May 

0.594 

0. 367 

0.032 

.lun . 

0.825 

0.427 

0.043 

.lul  . 

0.948 

0.563 

0.068 

Auj’ . 

0.607 

0.466 

0.052 

Sep. 

0.826 

0.505 

0.046 

0c  t . 

1.273 

0.494 

0.040 

Nov . 

1 . 192 

0.606 

0.044 

Dec . 

0.869 

0.554 

0.048 

Model 

ii 

m = 

.879  - 

. 050M  + 

. 006M* 

s = 

.282  - 

. 0 1 9M  l 

.00IM2 

KMS  Krrni  . 104 

Mode  1 

III 

. 002M  ' 

m * 

1.093  - 

. 2 1 5M 

+ .036M2 

s * 

.405  - 

. 07 6M 

+ .018M2  - 

. 00IM 

KMS  Error  * .101 

Mode  1 

IV 

m * 

.914  - 

. 06  5M  + 

.007M2  - . 

.179  Sin 

(n(M-3.5)/2) 

s * 

.272  + 

.02  1M  + 

.001M2  - . 

.037  Sin 

(ti  (M- 1.5)) 

KMS 

Krror  .082 

Mode  1 

V 

m * 

1.128  - 

. 2 30M 

+ .038M2  - 

. 002M  ’ + 

. 1 79  Sin  (n(M-l .5 

)/2) 

s * 

.364  - 

. 050M 

+ .014M2  - 

. 001M  - 

.032  Sin  (n(M-l  .5 

)) 

RMS  Krror  .079 


Station  14702  Bedford,  Mnssnclmsso t s 


Model  I 


m 

s 

RMS 

Jan. 

0.914 

0.486 

0.036 

Feb. 

0.902 

0.462 

0.024 

Mar. 

1.035 

0.521 

0.024 

Apr . 

0.959 

0.445 

0.050 

j 

May 

0.  787 

0.449 

0.051 

| 

Jun . 

0.739 

0.388 

0.031 

‘ 

Jul. 

0.861 

0.624 

0.023 

; 

Aug . 

0.866 

0.656 

0.046 

J 

Sep. 

0.788 

1.234 

0.055 

j 

Oct. 

0.993 

0.542 

0.052 

< 

Nov. 

1.386 

0.493 

0.051 

i 

Dec . 

0.970 

0.530 

0.050 

j 

Model  II 

m = 1.060 

- . 080M 

2 

+ . 007M 

i 

s = .324 

+ . 069M 

- .004M2 

RMS  Error  = .087 

Model  III 

« 

m = 1.003 

- . 037M 

- .001M2 

s = .755 

- . 264M 

+ .058M2  - 

. 003M 

3 

RMS  Error  = .085 

Model  IV 

m = 1.098 

- . 086M 

+ .007M2  + 

.100 

Sin  (if  (M-2.5)/2) 

s = .185 

+ .100M 

- .005M2  - 

.199 

Sin(i>  (M-3 . 5)  /3 . 5) 

RMS 

Error  = .075 

Model  V 

m = 1.186 

- . 1 51 M 

+ .019M2  - 

. 00  1 M 

5 + .112  Sin(ir(M-2.5)/2) 

s = .814 

- . 31 5M 

+ .068M2  - 

. 004M 

+ .101  Sln(ir  (M-l.  5) / 2 

.5) 

RMS  Error  = .076 
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Station  22521  Honolu lu , Hawa i i 
Model  I 


m 

s 

RMS 

Jan. 

0.947 

1 . 585 

0.048 

Feb. 

r>.  580 

0.662 

0.040 

Mar , 

0.444 

0.692 

0.068 

Apr. 

0.233 

0.349 

0.071 

May 

0.  189 

0.321 

0.088 

Jun . 

0.065 

0 . 099 

0.072 

Jul . 

0.148 

00 

o 

0.063 

Aug. 

0.122 

0.  163 

0.07  5 

Sep. 

0.205 

0.260 

0.044 

Oct . 

0.389 

0.496 

0.042 

Nov . 

0.558 

0.145 

0.056 

Dec . 

0.750 

1 .020 

0.046 

Model 

ii 

m = 

1.202 

- . 33 1M 

+ .025M2 

s = 

1.772 

- . 463M 

= .032M2 

RMS  Error  = 

.091 

Model 

III 

m = 

1.271 

- . 384M 

+ .035M2  - 

.OOIM* 

s = 

1.926 

- . 582M 

+ .054M2  - 

.00IM5 

RMS  Error  = 

.088 

Model 

IV 

m - 

1.195 

- . 330M 

+.025M2  - 

.025  Sin(u 

(M-l .5)) 

s = 

1.767 

- . 462M 

+.032M2  - 

.015  S i n ( it 

(M-l  . 5)) 

RMS 

Error 

Model 

V 

m = 

1.243 

- . 366M 

+.031M2 

- 

.022  Sin(n(M- 

1.5)) 

s - 

1.916 

- . 576M 

+.053M2  - 

.001M3  - 

.008  Sin  (it  (M- 

1.5)) 

RMS  Error  = 


'• 
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Station  26433  Ncjiana , Alaska 
Mode  1 I 


m 

s 

RMS 

Jan. 

0.  120 

0.  106 

0,039 

Feb . 

0.103 

0.1  14 

O 

o 

Mar. 

0.081 

0.046 

0.066 

Apr. 

0.075 

0.073 

0.046 

May 

0.  180 

0.117 

U.044 

Jim . 

0.428 

0.267 

0.026 

Jul  . 

0.415 

0.244 

0.048 

Auk- 

0.444 

0.314 

0.055 

Sep. 

0.272 

0.188 

0.043 

Oc  t . 

0.  124 

0.078 

0.032 

Nov . 

0. 103 

0,086 

0.023 

Dec. 

0.124 

0.114 

0.066 

Model  1 I 

m = - . 

105  + 

. I20M 

.009M2 

s = -. 

015  + 

. 06  1 M 

. 004 M 2 

RMS  Error  = .232 

Model  III 

m = .083  - . 02  5M  +.018M'  - .001M* 
s .103  - . 031M  +.013M2  - .001M3 

RMS  Error  = .220 

Mode  I IV 

tn  - -.077  + . 103M  - ,00/M2  - .146  Sin(n(M-2  )/3.5) 

s = -.076  + . 08 1 M - .006M2  - .089  Sin(ir (M-2 , 5) /3) 

RMS  Error  = .077 

Model  V 

m = -.096  + .072M  + .004M2  - .001M5  - .132  Sin  (it  (M-2 . 5) /3) 

s = -.009  + . 03 1M  + .004M2  - .00IM3  - .084  Sin (u (M-2 . 5) /3) 

RMS  Error  = .071 


>tat ion 

33123 

Tripol  i , Libya 

Mode  1 

J 

m 

s 

KMS 

Jan. 

0.434 

0.279 

0.053 

Feb. 

0.242 

0.475 

0.027 

Mar . 

0.167 

0.253 

0.034 

Apr. 

0 . 089 

0.214 

0.032 

May 

-0.008 

0.  123 

0.065 

Jun . 

-0.010 

0.028 

0.052 

Jul  . 

-0.076 

0 . 04  5 

0.016 

Aug. 

-0.034 

0.038 

0.01  5 

Sep , 

-0.394 

0.958 

0.040 

Oct . 

0.346 

0.387 

0.048 

Nov . 

0.679 

0.689 

0.035 

Dec . 

0.712 

0.465 

0.056 

Mode  1 

JJ 

m = 

. 777  - 

. 292M  4-  .024M2 

s = 

.468  - 

. 112M  + .OHM2 

KMS  Krror  = .135 

Model 

,002M3 

m = 

.474  - 

. 058M  - .019M2  + . 

s = 

. 723  - 

. 309M  + .047M2  - . 

, 002M  ! 

KMS  F.r ror  . 1 26 

Model 

JV 

m = 

.856  - 

. 304M  + .024M2  + 

. 146  S i 1 1 ( ti 

(M-3. 5)/ 3) 

s = 

. 243  - 

. 030M  + .005M2  - . 

.101  Sin(n(M-4. 5)/6) 

KMS 

Krror 

Model 

V 

m = 

.407 

- .014M  - .027M2  4- 

. 003M  ' - 

. 1 06  S in (n ( M-2 . 5) / 1 . 

5) 

s = 

1.012 

- . 627M  + . 116M2  - 

. 006M * 4 

.266  S i n ( u ( M- 1 .5) / 4) 

KMS  Krror  = .13  3 


I I-  I I 
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Station  41 108  Saigon , Vietnam 
Model  I 


m 

s 

RMS 

Jan . 

-0.435 

0.828 

0.055 

Feb. 

-0.081 

0.116 

C . 0 1 8 

Ma  r . 

-0.074 

0.457 

0.089 

Apr . 

0.234 

0.  592 

0.081 

May 

1.562 

2.184 

0.078 

Jun . 

2.201 

0.812 

C .081 

Jnl  . 

1.893 

0.787 

0.044 

Aug. 

2.004 

1.220 

0.077 

Sep. 

2.304 

0.681 

0.096 

Oct. 

2.073 

1.194 

0.041 

Nov . 

0.  761 

0.468 

0.052 

Doc . 

0.  345 

0.274 

0.015 

Mode  1 

j L 

m - 

-1.991 

F I.044M  - .069M2 

s = 

.053 

t . 322M  - .025M2 

RMS  Error  = .210 

Mode  1 

IN 

m = 

.579  - 

. 04 7M  + . 133M2  - . 

,010MJ 

s = 

.381  + 

. 069M  + . 022M2  - . 

.002M  ! 

RMS  Error  - .169 


Model 

IV 

m = 

2.202  + 1 

.07/M  - 

.069M2  - 

.546  Sin(ir(M-2.5)/2) 

s = 

.025  + 

. 326M  - 

. 025M-)  - 

. 102  Sin(n(M-2.5)/1.5) 

RMS  Error 

Modi*  1 

V 

m = 

-1.268  + 

. 38 5M  + 

.058M2  - 

. 00 7 M ^ - .421  S In  (it  (M-2 . 5 ) / 2 ) 

s = 

.235  + 

. 165M  + 

. 004M2  - 

.002M  1 - .232  Sin(ir  (M—  1 . 5) ) 

RMS  Error  = .154 
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v 





Station  45715  Shemya,  A I a ska 
Model  1 


m 

s 

RMS 

Jan. 

0.397 

0.213 

0.039 

Feb. 

0.101 

0.  197 

0.051 

Mar. 

0.304 

0.213 

0.042 

Apr . 

0.326 

0.292 

0.062 

May 

0.390 

0.233 

0.026 

Jun. 

0.325 

0.356 

0.045 

Ju  1 . 

0.530 

0.335 

0.042 

Aug. 

0.679 

0.355 

0.043 

Sep . 

0.588 

0.443 

0.038 

Oct . 

0.628 

0.340 

0.053 

Nov . 

0.542 

0.393 

0.030 

Dec. 

0.402 

0.296 

0.034 

Mode  1 

1 1 

m = 

. 215 

+ 

. 058M 

— 

.003M2 

s = 

. 120 

+ 

. 050M 

- 

.003M2 

RMS  F.rror  = .108 

Mode  1 

11  I 

m = 

. 584 

- 

. 227M 

+ 

.050M2  - 

. 003M * 

s = 

.235 

. 038M 

+ 

.014M2  - 

. 00 1M * 

RMS  I’.rror  = .059 

Model 

IV 

m = 

.299 

+ 

. 04  5M 

- 

.003M2  - 

. 145  S i n (tt (M—  1 . 5)  / 5) 

s = 

.118 

+ 

. 051M 

- 

.003M2  - 

.008  Sin(ir(M-l.  5)) 

RMS  F.rror 

Model 

_V 

m = 

.583 

— 

. 226M 

+ 

.050M2  - 

.003M3  + 

.001  S in (n (M-2 . 5) ) 

s * 

.232 

_ 

. 037M 

+ 

.01 3M2  - 

. 00 1 M 3 - 

. 002  S in(n(M- 1 . 5) ) 

RMS  Hr for  = 


11-13 


.062 


. 059 
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